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[i] Deep convective systems (DCSs) consist of intense convective cores (CC), large 
stratiform rain (SR) regions, and extensive nonprecipitating anvil clouds (AC). This 
study focuses on the evolution of these three components and the factors that affect 
system lifetime and AC production. An automated satellite tracking method is used in 
conjunction with a recently developed multisensor hybrid classification to analyze the 
evolution of DCS structure in a Lagrangian framework over the central United States. 

Composite analysis from 4221 tracked DCSs during two warm seasons (May-August, 

2010-2011) shows that maximum system size correlates with lifetime, and longer- lived 
DCSs have more extensive SR and AC. For short to medium systems (lifetimes <6 h), 
the lifetime is mainly attributed to the intensity of the initial convection. Systems that 
last longer than 6 h are associated with up to 50% higher midtropospheric relative 
humidity and up to 40% stronger middle to upper tropospheric wind shear. Such 
environments allow continuous growth of detrained hydrometeors by deposition, 
supporting further development of the SR and AC region, as indicated by the increased 
staggered timing between stratiform clouds and peak convective intensity, thus 
prolonging the system lifetime beyond 6 h. Regression analysis shows that the areal 
coverage of thick AC is strongly correlated with the size of CC, updraft strength, and SR 
area. Ambient upper tropospheric wind speed and wind shear also play an important role for 
convective AC production, where for systems with large AC (radius >120 km) they are 24% 
and 20% higher, respectively, than those with small AC (radius = 20 km). 

Citation: Feng, Z., X. Dong, B. Xi, S. A. McFarlane, A. Kennedy, B. Lin, and P. Minnis (2012), Life cycle of midlatitude deep 
convective systems in a Lagrangian framework, J. Geophvs. Res., 117, D23201, doi:10.1029/2012JD018362. 


1. Introduction 

[ 2 ] A deep convective system (DCS) is a contiguous high 
cloud shield that consists of active deep convective cores 
and citrus anvil canopy. The convective cores produce 
heavy precipitation and provide moisture to the upper tro- 
posphere to maintain the DCS stratiform precipitation and 
the nonprecipitating anvil cloud [Lin et al, 2007]. The anvil 
cloud shields dominate the DCS radiative budget due to their 
large spatial coverage. 

[ 3 ] Because they cover large areas, satellite observations 
are needed to examine the radiative impact of DCSs and to 
study the life cycle of convective clusters [. Machado and 
Rossow, 1993; Del Genio and Kovari, 2002; Feng et al.. 
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2011]. Williams and Houze [1987] developed an objective 
tracking technique using geostationary satellite data to study 
the cloud clusters over the western Pacific maritime conti- 
nent during the Winter Monsoon Experiment. Pope et al. 
[2008] used multiyear Japanese Meteorological Agency 
Geostationary Meteorological Satellite-5 Infrared (IR) data 
to track and study the evolution of Mesoscale Convective 
Systems (MCSs) over northern Australia and surrounding 
oceans. They found that the minimum brightness tempera- 
tures of MCSs decreased slowly as their sizes reached 
maximum values during their mature stage, suggesting that 
the areal growth during that stage results from anvil forma- 
tion, although they could not determine the amount of cirrus 
anvils by using IR data alone. Futyan and Del Genio [2007] 
used geostationary Meteosat-8 satellite data to identify and 
track organized convective systems over tropical African 
and Atlantic regions, and subsequently matched these sat- 
ellite data with overpasses from the Tropical Rainfall Mea- 
surement Mission (TRMM) satellite and composited vertical 
structure based upon their related life cycle stages. They 
found that African storms are convectively active with fre- 
quent lightning during their developing stages and evolve 
into more stratiform structures as they dissipate. In general, 
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intense storms have larger cloud sizes, greater thicknesses, 
and longer durations than weak storms. They hypothesized 
that convective intensity might be a useful statistical pre- 
dictor for use in generation circulation models (GCMs) to 
simulate the eventual sizes of DCSs without resolving 
mesoscale dynamics. 

[ 4 ] In the midlatitude, particularly in the United States, 
DCSs occurring in the warm season (May-August) are often 
associated with a wide variety of severe weather events and 
significant precipitation. Maddox [1980] first noted a par- 
ticular type of long-lasting MCS, termed Mesoscale Con- 
vective Complex (MCC) that frequently occurs over the 
central United States, grows to tremendous areal extent and 
often persists over 12 h. Composite analysis of satellite 1R 
and surface precipitation data by McAnelly and Cotton 
[1989] showed that average rainfall intensity peaks early in 
the satellite IR-defined life cycle, similar to the results 
shown by Kane et al. [1987], Jirak et al. [2003] compared 
several types of MCSs using satellite and radar data during 
the warm seasons and found that elongated MCSs are most 
common, and are most likely associated with severe weather 
and the most intense precipitation. In additional to their 
impact on society through severe weather, these systems also 
modulate the climate system through their impacts on the 
water cycle and cloud radiative forcing and feedback. 
Despite their importance, most GCMs do not correctly 
simulate the diurnal cycle of warm season continental pre- 
cipitation [Dai, 2006], in part due to the absence of meso- 
scale organization that generates extensive stratiform rain, 
preventing the lifetime of DCSs from extending beyond the 
decay time of individual convective events [ Del Genio et al, 
2012]. This error could have significant impact on the cloud 
radiative feedback associated with DCSs in GCMs. Feng 
et al. [2011] found that stratiform rain and anvil clouds 
from DCSs in the southern great plains of the United States 
have peak areal coverage from late evening to early morning, 
and that anvil clouds are major contributors to the summer- 
time TOA radiative budget. Errors in the timing and extent of 
DCS stratiform rain and anvil clouds may result in incorrect 
diurnal forcing due to erroneous cloud radiative feedbacks. 

[5] This study focuses on the DCS cloud life cycle in the 
central United States, particularly on the evolution from 
convective to stratiform anvil cloud structures, and the fac- 
tors controlling convective anvil production. Previous stud- 
ies of DCSs over various tropical land locations found many 
similarities in terms of cloud life cycles [Machado et al, 
1998; Machado and Laurent, 2004; Kondo et al, 2006; 
Futyan and Del Genio, 2007; Pope et al, 2008], If these 
characteristics hold true for midlatitudes as well, such 
“universal” behavior would provide useful insights and 
implications for parameterization development in future 
GCMs. The large spatial coverage of the Next-Generation 
Radar (NEXRAD) network and routinely available Geosta- 
tionary Operational Environmental Satellite (GOES) data 
over the continental U.S. are valuable observational resour- 
ces for studying the DCS life cycle in the midlatitudes. 
By implementing a satellite cloud tracking method along 
with the hybrid classification developed by Feng et al. [2011], 
the evolution of the major DCS components is analyzed here 
in a Lagrangian framework to study the morphology of mid- 
latitude DCS clouds. Section 2 describes the data sets and 
methodology used in this study; section 3 discusses the 


evolution of DCS structures. Section 4 investigates factors that 
affect anvil production, and finally, a summary and concluding 
remarks are given in section 5. 

2. Data and Methodology 

[6] The data sets used in this study are similar to those in 
Feng et al., [201 1]: both NEXRAD radar and GOES satellite 
observations over the central United States (25-40°N, 80- 
105°W) during two warm seasons (May-August) in 2010 
and 2011. However, the focus of this study is on the DCS 
life cycle, while the previous one was on the TOA radiation 
budget. 

[7] The NEXRAD data were obtained from the National 
Mosaic and Multisensor Quantitative Precipitation Estimate 
(QPE) system (NMQ) [Zhang et al., 201 1]. Reflectivity data 
from all NEXRAD radars in the contiguous United States 
were quality controlled and objectively projected onto a 
unified three-dimensional Cartesian grid. Instantaneous rain 
rate data provided in the QPE product are derived using Z-R 
relationships for convective, stratiform and tropical rain 
types. The NMQ three-dimensional Mosaic reflectivity data 
set has a fixed 0.01° (~1 km) horizontal resolution, a 0.5 km 
vertical resolution up to 18 km above mean sea level, and 
5 min temporal resolution. The GOES cloud products are 
retrieved using the algorithms developed for the NASA 
Clouds and Earth’s Radiant Energy System (CERES) 
project [Minnis et al., 2008, 2011], The GOES cloud 
properties used in this study are 10.8 //m channel bright- 
ness temperature (T IR ), and cloud top and cloud base 
heights with a pixel resolution of about 4.5 x 10 km and 
temporal resolution of 30 min. 

[8] The three-dimensional temperature, moisture, and 
wind fields used in this study were obtained from the 
National Centers for Environmental Prediction (NCEP) 
North America Regional Reanalysis (NARR) [Mesinger 
et al., 2006]. The NCEP NARR data set has 3 h temporal, 
32 km horizontal, and 45-layer vertical resolutions over the 
North American domain. NARR uses the operational NCEP 
Eta Model and its three-dimensional variational data assim- 
ilation technique to improve on NCEP reanalysis by accu- 
rately assimilating precipitation [Kennedy et al., 2011]. The 
use of the NARR data set is intended to investigate large- 
scale circulation and moisture variations associated with 
DCS development, rather than resolving convective-scale 
features (e.g., updraft velocity) due to its spatial and tem- 
poral resolution. 

[ 9 ] The multisensor hybrid classification technique 
developed by Feng et al. [2011] is applied to the matched 
NEXRAD and GOES data sets. This technique jointly ana- 
lyzes the collocated three-dimensional NEXRAD reflectivity 
structure and GOES T IR fields to objectively identify DCSs 
and to subsequently classify convective core (CC), strati- 
form rain (SR) and anvil cloud (AC) components from 
DCSs. Compared to previous studies that use only satellite 
T ir thresholds to define DCS and separate convective core 
and anvil clouds, the hybrid classification method avoids 
misidentification of nonconvective cold cirrus clouds as 
DCS, and explicitly separates SR, thick and thin nonraining 
ACs from CC. This is particularly suited for investigating 
their life cycles because the dynamical and microphysical 
properties are distinctly different between these regions. The 
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1 km resolution NEXRAD radar data are convolved into the 
GOES pixels as described by Feng et al. [2011], 

[ 10 ] To obtain the environment where DCSs form, the 
3 -hourly NARR data set is used at the time when it is 
available and collocated with DCSs identified by the 
hybrid classification. For example, if a DCS initiates at 
00 UTC and dissipates at 08 UTC within the study 
domain, the environmental condition for this DCS is then 
obtained from NARR at 00, 03 and 06 UTC. While NARR is 
not available at each hourly time step of all the DCSs, by 
compositing large data samples the overall environmental 
condition associated with DCSs of a given lifetime or AC 
area can be obtained. The time-matched and collocated 
NARR wind and moisture profiles are averaged within indi- 
vidual DCSs and saved for further analysis. 

[n] An automated satellite tracking method is developed 
to study the life cycles of DCSs by tracking the evolution of 
the systems and subsequently incorporating with the 
corresponding classification of CC, SR, and AC regions 
from the hybrid classification. This tracking method, similar 
to that described by Williams and Houze [1987] and Futyan 
and Del Genio [2007], is based on coincidence overlapping 
of successive hourly GOES T IR data. Each contiguous 
cloudy region with T IR <235 K is used to define a cold cloud 
system for tracking. Only cold cloud systems with a mini- 
mum area of 400 km 2 (~8 GOES pixels) are tracked by the 
algorithm. A secondary T 1R threshold (215 K) is used to 
identify areas of strongly depressed IR temperature, labeled 
as cold cores. Both cold cores and cold clouds are used for 
tracking the convective systems. If, within two successive 
T ir images, a cloud system overlaps itself by more than 50% 
of its area, this cloud system is considered as the same sys- 
tem and its track continues. Tracking is performed iteratively 
on successive hourly T IR images. Sensitivity tests show that 
changing the cold cloud T IR and area overlap thresholds do 
not affect any conclusions. 

[ 12 ] The overlap comparison is performed either going 
forward or backward in time for all cloud systems that have 
overlaps. If a cloud system meets the overlap requirement in 
either tracking direction, it is considered as the same system. 
When two or more systems merge together, the largest 
overlapped system is assumed to continue and the smaller 
overlapped systems are terminated. Similarly, when a cloud 
system splits into several smaller systems, the largest over- 
lapped fragment continues carrying the characteristics of the 
original system and the smaller overlapped fragments are 
labeled as newly formed clouds. The merging and splitting 
processes are recorded for each time step, thus the systems 
that experience merging or splitting during their life cycles 
can be identified in further analysis. 

[ 13 ] To track the DCSs through their life cycles and to 
incorporate the large areal coverage of AC associated with 
DCSs, the tracked cold clouds are subsequently expanded 
outward to T IR up to 270 K, as defined by the hybrid clas- 
sification [Feng et al., 2011]. By matching the tracked cold 
clouds with the hybrid classification, the life cycle of pre- 
cipitation and cloud structures associated with each DCS are 
obtained for the following analysis. 

[ 14 ] Figure 1 illustrates an example of the life cycle of a 
DCS tracked by the automated tracking method on 14-15 
June 2010, with the corresponding classification of CC, SR, 
and AC regions. Note that AC represents all three anvil 


subcategories (i.e., transitional, thick, thin anvils) as defined 
by both NEXRAD and GOES data in the hybrid classifica- 
tion [Feng et al., 201 1], and characterizes the majority of the 
anvil clouds associated with the DCS. The system was ini- 
tially detected at 1715 UTC and grew in size continuously 
for the following 10 h. At 2215 UTC, the system began 
separating into two parts and eventually at 0015 UTC of 
15 June 2010, the relatively smaller system over eastern 
Oklahoma split from the main system and started a new 
track system at 0 1 1 5 UTC (not part of this system and hence 
not shown). The mature stage of the system (02 1 5-061 5 UTC, 
15 June 2010) represented a squall line structure of leading 
CC (red color) and trailing SR (green color), with extensive 
area of associated peripheral AC (blue color). The system 
started dissipating at around 0715 UTC and the tracking was 
eventually terminated at 1315 UTC of 15 June 2010. This 
particular system developed and matured as an MCS while 
slowly propagating eastward between the boundaries of Texas 
and Oklahoma during its life cycle, and was well captured by 
the automated tracking method. 

[ 15 ] Radar reflectivity has been widely used to indicate the 
strength of convective updrafts [Houze et al., 1990; Steiner 
et al., 1995; Zipser and Lutz, 1994], For each tracked sys- 
tem, several convective intensity (Cl) indices are calculated 
from the NMQ radar data set, namely the average instanta- 
neous precipitation rates of convective core ( PR C c ) and 
stratiform rain regions ( PRsr ), the maximum height of the 45 
dBZ contour in CC (H 45dBZ , proxy for strong updraft as 
shown by Zipser et al. [2006]), and equivalent radius of the 
area with reflectivity greater than 30 dBZ at 8 km height 
C^ccs-km, proxy for the size of intense convective cores, 
same as that used by Rickenbach et al. [2008]). 

[ 16 ] Figure 2 shows an example of the evolution of a 
tracked DCS and its associated cloud top temperature, CC, 
SR, and AC properties obtained from the hybrid classifica- 
tion during its life cycle. Three stages (developing, mature 
and dissipating) are defined using the evolution of the DCS 
equivalent radius and IR temperature following the life cycle 
stage definition of Futyan and Del Genio [2007]. When a 
system is detected and tracked, the system is considered as 
growing vertically and defined as “developing” before it 
reaches the minimum IR temperature (Figure 2b). The sub- 
sequent period is considered as the “mature” stage when the 
system reaches the minimum IR temperature and continues 
to grow in size until reaching its maximum (Figure 2a). The 
“dissipating” stage is defined after the mature stage when the 
system decreases in size. 

[ 17 ] As shown in Figure 2c, CC area grew significantly 
and reached its maximum size (corresponding to the mini- 
mum IR temperature in Figure 2b) during the developing 
stage and then decreased slowly at the following two stages. 
SR and AC area coverage also increased rapidly during the 
developing stage, but they kept a more gradual expansion 
and reached their respective maximum sizes at the end of the 
mature stage, and then both decreased sharply in the dissi- 
pating stage. As a result of detrained ice and/or mixed phase 
hydrometeors from the CC region, the variations of the SR 
and AC area showed similar trends, and their evolutions 
were closely tied to each other, suggesting that SR and AC 
are a contiguous cirrus anvil cloud shield and only split by 
the analyzing technique for rain and no rain. While the SR 
region has precipitation detected by radars, hydrometeors in 
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Figure 1. The life cycle of a convective system that is tracked by the automated tracking algorithm in the 
study domain. Time increases from the top left to the bottom right, and each image represents an hour. The 
colors represent regions given by the hybrid classification. 


the AC region mostly remain aloft. The several hour lag 
between peak sizes of stratiform anvil region (both SR and 
AC) and CC in Figure 2c indicates the gradual transition 
from a mostly convective to a more stratiform structure for 
this particular system. Similar to the life cycle of CC size, 
both PRcc an d H 45 dBZ peaked at the early life stage, and then 
gradually decreased during the mature and dissipating stage 
(Figure 2d). This example suggests that the Cls examined 
here peak between developing and mature stages, about at the 
same time as the maximum CC size but 5-6 h ahead of the 
maximum SR and AC sizes. Following the same analysis 
method, the evolution of these properties from all tracked 
systems are statistically evaluated in the next section. 

3. Evolution of DCS Properties 

[is] A total of 10599 DCSs were tracked during the two 
warm seasons (May- August) in 2010 and 2011, resulting in 
an average frequency of 43 DCSs per day over the study 
domain. The average DCS cloud fraction accounts for 25% 
of total clouds observed by GOES satellite (8% of total 


clouds are from precipitating convective and stratiform rain 
portions of the DCS, and 17% of total clouds are from anvil 
clouds), or 47% of all high clouds, consistent with our pre- 
vious study [Feng et al., 2011]. In comparison, using 30 min 
GOES imagery to track convective systems, Bedka and 
Minnis [2010] found 877 DCSs in 22 days (~39 per day) 
during the 2007 Tropical Composition, Clouds and Climate 
Coupling Experiment over Central America. 

[ 19 ] To obtain robust statistics that accurately represent 
characteristics of DCSs throughout their life cycles, systems 
with centroids located within 1° latitude/longitude of the 
boundaries of the study domain during any time in their life 
cycle are not used. Further, if the maximum radar coverage 
area over the entire life cycle of a tracked DCS is less than 
70% of its size, the system is possibly outside the radar 
coverage area in the ocean or near the Rocky Mountain at 
the western edge of the domain. Precipitation and anvil 
structures of these systems are ambiguous and thus excluded 
from the analysis. With these criteria, the number of systems 
is reduced to 4221, or ~40% of total samples, and their 
PDFs and CDFs are shown in Figure 3. 
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Figure 2. Example of the evolution from a tracked deep convective system (DCS) showing (a) system 
radius, (b) minimum IR temperature, (c) radii of convective core (CC), stratiform rain (SR), and anvil 
clouds (AC) calculated from the hybrid classification, and (d) CC rain rate (PRccX 45 dBZ echo top height 
(H 45 dBz)- The radius is defined as the equivalent radius of a circle (R), such that the area is given by 
A = i xR 2 . Lines are forth-order polynomial fits, and symbols are observed values. Numbers on top of 
the panels give the life cycle stage: 1, developing; 2, mature; 3, dissipating. 


[ 20 ] As expected, system lifetimes decrease exponentially 
with increasing lifetime (Figure 3 a) with a dominant popu- 
lation of short-lived systems (lifetime <3 h, 49%), and fewer 
than 1% of the systems having a lifetime longer than 14 h. 
This distribution of lifetimes is comparable to that found by 
Bedka and Minnis [2010] for a domain centered on Central 
America. The longest-lasting systems (usually MCSs) con- 
tribute significantly to warm season precipitation in the 
central U.S. and thus are crucial to the hydrological cycle 
[' Carbone et al., 2002]. From the tracked DCS database, the 
longest-lived systems (>14 h) contributed ~27% rainfall and 
1 7% cloud cover, however, these systems were not included 
in the statistical analysis due to limited data samples. 

[ 21 ] Among 4221 analyzed DCSs, ^90% have an equiv- 
alent radius less than 120 km and a minimum cloud top T IR 
lower than 220 K (Figures 3b and 3c). Maximum convective 
core radius ( Rcc )> stratiform rain radius (Rsr) and anvil 
cloud radius (Rac) (Figures 3d— 3f) are obtained from the 
hybrid classification. Although the modal values of R cc and 
R S r are about the same (~20 km), larger SR occur more 
frequent than CC (e.g., 90th percentile of Rsr is at 60 km 
while it is at 35 km for R cc ). The average area covered by 
rain from CC and SR during the study period are 28% and 
72%, respectively, similar to the average tropics-wide SR 
area fraction (73%) documented by Schumacher and Houze 
[2003], Compared to the CC and SR coverage, the AC 
dominates the DCS area and its distribution (Figure 3f) is 
similar to that of the DCSs (Figure 3b). At the 90th per- 
centile, the area covered by AC is about three times that of 


SR and almost an order of magnitude that of CC. This result 
is consistent with Feng et al. [2011] over the southern great 
plains, suggesting that this attribute is representative of the 
DCSs in the southern and central U.S. region. For the DCSs 
with reflectivity greater than 30 dBZ at 8 km MSL height, 
most of their radii (7?cc8-km) are 10 km, with a rapid decrease 
to ~20 km at the 90th percentile (Figure 3g) due to few 
systems with broad areas of intense updrafts supporting 
large particles at this height. The modal value for maximum 
45 dBZ echo top height {H 45m7 ) is around 5-6 km, with the 
90th percentile at a 9 km height (Figure 3h). The distribu- 
tions of average CC precipitation rates ( PRcc > Figure 3i) are 
similar to those in the SR regions ( PRsr , not shown), except 
the magnitude of PRcc is an order of magnitude that of 
PRsr- The 90th percentile of PRcc is about 12 mm h -1 and 
1.4 mm h -1 for PRsr. 

[ 22 ] The relationship between total system lifetime and 
maximum equivalent radius for the tracked DCSs is illus- 
trated in Figure 4. Tracked DCSs that do not experience 
merging or splitting during their lifetimes are shown sepa- 
rately in blue colors. This analysis found that only ~4% of 
the tracked DCSs split into multiple systems. The maximum 
DCS size attained during DCS life cycles correlates with 
system lifetime. However, the variance for longer-lived 
systems increases significantly, similar to the results of 
Machado et al. [1998] and Pope et al. [2008]. This is pos- 
sibly due to limited samples observed for large and long- 
lived systems, resulting in larger statistical errors for the 
averages. When the correlation is calculated based on 
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Figure 3. PDF (bars) and CDF (black lines) of (a) system lifetime, (b) maximum size for each tracked 
DCS, (c) minimum IR temperature, (d) maximum convective core radius ( Rcc X ( e ) maximum stratiform 
rain radius ( R sr ), (f) maximum anvil cloud radius ( Rac X ( S ) radius of reflectivity greater than 30 dBZ area 
at 8 km height (Rccs-kmX (h) 45 dBZ echo top height (// 4 5 dBz)> and (i) convective precipitation rate 
(PRcc). 


lifetime <10 h (red dash line), the average radii for longer 
systems are still within the range of linear fit lines from the 
25th and 75th quartiles (shaded area). The maximum sizes of 
nonmerged systems also correlate with lifetimes up to a 
period of 10 h. The number of systems decreases dramati- 
cally for nonmerged systems at lifetimes longer than 10 h, 
suggesting that long-lasting DCSs usually experience 
merging with other systems. These results are consistent 
with those of Jirak et al. [2003] in that more than 70% of 
their MCS samples evolved from the merger of multiple 
convective clusters. Comparisons between merged and 
nonmerged systems in this study (not shown) suggest that 
the major differences are system lifetime and maximum size 


attained, while other Convective Intensity (Cl) indices are 
similar. 

[ 23 ] To examine the evolution of CC, SR and AC of the 
tracked DCSs, the equivalent radii of these three components 
are composited and illustrated in Figure 5 based on their 
lifetimes following the method used by Pope et al. [2008]. It 
is apparent that the maximum sizes of the three major DCS 
components attained during their life cycles also correlate 
with system lifetime, i.e., longer-lived DCSs are consistently 
larger in size with more extensive stratiform rain and anvil 
cloud shields. This correspondence between system lifetime 
and maximum size is consistent with findings of Chen and 
Houze [1997] and Pope et al. [2008], suggesting systematic 
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Figure 4. The maximum DCS equivalent radius as a func- 
tion of system lifetime. Red colors are all systems including 
merged systems, and blue colors denote nonmerged systems. 
The box covers the 25th-75th percentiles, and the whisker 
covers the 10th-90th percentiles. The square indicates the 
mean value, and the cross denotes the median value. The 
dashed lines are linear function fitting using mean values 
for all systems up to 10 h (red) and for nonmerged systems 
up to 7 h (blue), with corresponding correlation coefficients 
shown in parentheses. The shaded areas are the same linear 
fits except using the 25th and 75th percentile values, with 
dark gray from all systems and light gray from nonmerged 
systems. Solid circles and diamonds are the number of sys- 
tems for each lifetime bin. 


similarities between tropical and midlatitude DCSs. Com- 
paring the life cycles of Rcc (Figures 5a and 5d) with 
those from R sr and R AC (Figures 5b, 5c, 5e, and 5f), it 
appears that the size of CC peaks always lead those of SR 
and AC by a few hours for all lifetimes and system sizes. 
Since Rsr and R A( ; are highly correlated (correlations 
>0.95 for all systems with the range of lifetimes shown in 
Figure 5), R AC is used to represent the entire mesoscale 
stratiform cloud deck (SR + AC) of these systems. 

[ 24 ] To understand the phase relationship between con- 
vective intensity and stratiform cloud coverage, Figure 6 
shows the lag between peak Cl indices and maximum R AC . 
It is clear that this lag increases linearly with system lifetime, 
regardless of which Cl index is used although the increase 
calculated from R cc is ^40% slower than those from the 
other three CIs. The radar-reflectivity-texture-based CC/SR 
partitioning [Steiner et al., 1995] used in the hybrid classi- 
fication tends to produce excessive CC and low SR area in 
enhanced SR rain [Lang et al., 2003], thus delaying the peak 
in R cc . Nevertheless, all of the CIs indicate that the 
increasing lag with system lifetime is a robust signal. The lag 
increases from ~ 1 h for 4 h long systems to more than 3 h for 
14 h long systems. 

[ 25 ] Occurrence of the peak convective intensity and 
stratiform cloud region (SR and AC) portions of DCSs at 
different times has important implications for climate 
model simulations. In the life cycle of DCSs, precipitation 
evolves from predominately convective in the early stage 
to more stratiform in the mature to decaying stage. The 
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Figure 5. Averaged system evolution for the three major components of the tracked DCSs. (left) All sys- 
tems, including merged systems, and (right) systems that do not merge or split, (a, d) Convective core 
(CC), (b, e) stratiform rain (SR), and (c, f) anvil cloud (AC). Dotted lines are hourly mean observed values, 
and solid lines are third-order polynomial fits. 
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Figure 6. Composite lag between maximum stratiform rain 
(SR) or anvil cloud (AC) area and peak convective intensity 
indices as a function of system lifetime. 

accompanying diabatic heating structures shift from the 
first baroclinic mode, which has heating throughout the 
troposphere, to the second baroclinic mode, which has 
heating in the upper troposphere and cooling below 
[Houze, 1989]. This change in the heating structure gen- 
erates mesoscale circulations in the upper troposphere that 
helps maintain the systems, particularly for the stratifonn 
area (SR and AC), thus, prolonging the lifetime of sys- 
tems. Note that microphysics also plays an important role. 
For example, incorrect hydrometeor size distributions in 
either CC or SR regions, or poor particle fall velocity 
assumptions in model simulations will result in failure to 
reproduce the observed structures in both CC and SR 
region [Varble et al., 2011], and the time delay charac- 
teristic between the convective and stratifonn cloud 
region. 

[ 26 ] In a comparison of Australian monsoon season MCSs 
simulated by a cloud resolving model (CRM) during the 
active and break periods of the Tropical Warm Pool-Inter- 
national Cloud Experiment (TWP-ICE), Del Genio et al. 
[2012] found ^4 h difference in lags between the SR 
region and convection during the two periods. They 
hypothesized that in a drier and less convectively sustainable 
environment, SR area is restricted to the amount of hydro- 
meteor detrainment from CC, and the detrained ice sublimes 
with no development of mesoscale circulations. In contrast, 
in a moist environment, upper tropospheric relative humidity 
(RH) is almost water saturated, and detrainment of buoyant 
air initiates a mesoscale updraft in the upper troposphere. 
High RH provides an environment favorable for sediment 
ice to grow further by deposition, releasing latent heat and 
eventually generating mesoscale updraft that maintains SR 
for hours. The much longer lag time represents the time it 
takes for mesoscale updraft to develop and maintain the SR 
region. 

[ 27 ] According to this hypothesis, the increased lag 
between maximum SR area and peak Cl can be explained as 
follows. In an environment favorable (or “sustainable”) for 
convective development (e.g., deeper moisture, higher 
CAPE, stronger wind shear, and/or presence of a low-level 
jet as suggested by Coniglio et al. [2010]), longer-lived 
DCSs/MCSs are more likely to form than shorter-lived sys- 
tems. Once convection develops, ambient upper level wind 


spreads out the initial detrained ice from the CC, and the 
humid upper tropospheric environment allows development 
of mesoscale updrafts through deposition growth of sedi- 
ment ice that releases latent heat. This process takes several 
hours until the mesoscale stratiform cloud region is fully 
developed. As a result, the lag between maximum SR and 
AC area and peak Cl for long-lived systems is considerably 
larger than that for the short-lived ones. The relationship 
between system lifetime, convective organization and the 
environment will be further investigated in section 4. 

[ 28 ] Note that the 1-3 h lags in Figure 6 represent the 
composite mean of all tracked systems. A 1-2 h lag between 
maximum convective rainfall and system cloud cover was 
reported by Rickenbach et al. [2008] for several isolated 
convective cells in south Florida. The dry and suppressed 
environment in their study and lack of mesoscale organiza- 
tion represents the “unsustainable” environment and resem- 
bles the low end of the lag. Individual DCSs in this study 
show much larger spread of lags. The calculated standard 
deviation of lags has a range of 0-7 h, with the upper end 
comparable to the TWP-ICE moist case shown by Del Genio 
et al. [2012], Nevertheless, the composite lags show the 
increased staggered timing between convective and strati- 
fonn latent heating, and radiative heating from anvil clouds. 
Latent heating profiles associated with convective clouds are 
positive throughout the troposphere due to active conden- 
sation of droplets, while stratifonn profiles have heating 
above the melting level due to condensation/deposition and 
cooling below due to evaporation of rain [Schumacher et al., 
2004]. Radiative heating profiles from optically thick AC are 
positive at AC base (between middle and upper troposphere) 
and negative at AC top, while optically thin cirrus have 
radiative heating throughout the cloud [Ackerman et al., 
1988; Garrett et al., 2005; Jensen and Del Genio, 2003; 
Powell et al., 2012], As a result of DCSs transitioning from 
convective to stratifonn structure, the positive latent heating 
shifts to the upper troposphere, and is further enhanced by 
the radiative heating from AC due to the concunent peak of 
the SR and AC regions. The increased staggered timing with 
overall system lifetime suggests that changes in the overall 
diabatic heating structure could help maintain the systems 
through mesoscale ascent, prolonging the system lifetime by 
sustaining SR and AC development. 

[ 29 ] The lack of mesoscale organization in most current 
generation GCMs prevents them from even forming a SR 
region [Del Genio et al., 2012], let alone correctly simulat- 
ing the staggered timing. CRMs have the ability to produce 
mesoscale features associated with DCSs through explicit 
simulation of cloud microphysics, although they have not 
yet produced realistic radar observed stratifonn structures 
[Varble et al., 2011], The evolution of the CC/SR/AC 
regions and their staggered timing, found in this study, 
should be compared with more CRM simulated cases to 
diagnose the mechanisms responsible for transition from CC 
to SR and AC, and provide further guidance for parameter- 
izing these processes in GCMs. 

4. Factors Affecting Anvil Production 

[ 30 ] For GCMs to predict how much stratiform rain (SR) 
and anvil clouds (AC) are produced from deep convection, 
it is important to determine the factors affecting their 
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Figure 7. Frequency distribution of the maximum anvil cloud radius as a function of various convective 
intensity (Cl) indices during developing and mature stages, (a) CC radius, (b) SR radius, (c) 45 dBZ echo 
top height, (d) radius of 30 dBZ area at 8 km height, (e) CC rain rate, and (f) SR rain rate. Black dots rep- 
resent median AC values for each Cl bin, the dashed line is the linear fit line from all raw samples, and r is 
the correlation coefficient. Results are computed using systems that do not experience merging or splitting 
during their life cycles. 


production so that these relationships can be correctly 
parameterized in the models. We explore the relationship 
between the size of AC and various Convective Intensity 
(Cl) indices, as well as the environment, with the intent of 
determining which factors are most relevant to the con- 
vective AC production over the study domain. To avoid 
influence from processes not related to the internal 
dynamics, the results in this section only use systems that 
did not experience any merging or splitting during their 
entire life cycles. The number of systems is reduced to 


3450 (81.7% of the original sample number), and most of 
these systems have lifetimes shorter than 10 h (Figure 4). 

[ 31 ] Figure 7 shows the frequency distributions of AC 
radius as functions of a variety of DCS properties. AC radius 
(Rac) and DCS properties are reported as the maximum 
values obtained during the developing and mature stages of 
the systems (see Figure 2). Similar to the results found in 
Figure 5, R A c correlates well with CC radius (Rcc) and SR 
radius (R,s/t) (Figures 7a and 7b), with correlation coeffi- 
cients of 0.50 and 0.64, respectively. The strong correlation 
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between AC and SR is expected because SR and AC are 
essentially a contiguous stratiform cloud deck from DCSs. 
Similar correlations between convective area and anvil area 
were found by Li and Schumacher [2011] in several tropical 
regions using TRMM precipitation radar data. However, 
anvil area was only weakly correlated with stratiform rain 
area in their study, possibly due to the low sensitivity of 
TRMM precipitation radar to nonprecipitating anvil clouds 
or to the narrow TRMM radar swath, which may not have 
completely observed some very large DCSs. 

[ 32 ] The correlation between R AC and several other Cl 
indices (Figures 7c-7f), namely the 45 dBZ echo top height 
(# 45 dBzX the radius of 30 dBZ area at 8 km (Rcc-8kmX CC 
precipitation rate ( PRcc ) and SR precipitation rate ( PRsr ), 
are much lower than those of Rcc and Rsr- The entire anvil 
cloud shield associated with a mature DCS (i.e., R AC ) con- 
sists of two parts: thick AC consisting of large hydrometeor 
particles that is still attached to CC or SR areas, and thin AC 
consisting of cirrus that is a remnant of the SR or thick AC 
region [. Houze , 1993; Li and Schumacher, 2011]. In the 
hybrid classification, the former is classified by NEXRAD 
radar, and the latter is identified by the GOES data. On 
average, thick anvils are 27% of the total anvil cloud fraction 
during the study period. One possible reason for the low 
correlation between total AC size and Cl is that thin AC, 
usually located furthest away from the active precipitating 
area, is less directly impacted by the parent convection and 
more likely influenced by the environment such as wind 
shear, moisture and radiative forcing from the atmosphere 
and surface underneath [ Garrett et al., 2005; Rickenbach 
et al., 2008]. 

[ 33 ] To investigate this possibility, each of the Cl indices 
was correlated with the radius of thick AC (not shown). The 
resulting correlation coefficients with H 45 dBZ , R C c- 8 km- and 
PRcc increased to 0.40, 0.59 and 0.35, respectively. The 
multivariate correlation coefficient between all CIs (shown 
in Figure 7) and the thick AC also increased to 0.79, com- 
pared to 0.69 for the total AC. Therefore, the most important 
factors that control the areal extent of thick AC are: the size 
of CC (Rcc, ^cc- 8 km) and the strength of convective updraft 
(proxies by // 45 dBZ> PRcc)- De/ Genio et al. [2012] noted 
that cumulus mass flux, defined as convective area times the 
mean updraft speed, is the best predictor for the SR area in 
their model simulation. The result in this paper is qualita- 
tively consistent with that suggested by Del Genio et al. 
[ 2012 ], 

[ 34 ] To investigate the relationship between the environ- 
ment and AC production, vertical profiles of atmospheric 
state variables from the NARR reanalyses were matched 
with the tracked DCSs. The results from NARR were cal- 
culated using data between May and August 2010 because 
the 201 1 NARR data set has not been released at the time of 
this paper. Figure 8 shows the composite profiles of zonal 
and meridional wind, relative humidity (RH), and vertical 
motion as a function of the maximum AC area attained 
during the developing and mature stages of the tracked 
DCSs. It is evident from the zonal wind composite 
(Figure 8a) that DCSs with larger AC area are associated 
with stronger upper level wind between 150 and 200 mb and 
deep layer wind shear. For example, the average wind speed 
between 150 and 200 mb for systems with R AC >120 km is 


24% higher than those with R AC = 20 km. The upper tro- 
pospheric bulk wind shear (wind speed difference between 
150 mb and 600 mb) increases ~20% monotonically from 
13.2 m s -1 (for R AC = 20 km) to 15.9 m s -1 (for R AC >120 
km). The evidence of stronger upper level wind speed and 
deep layer wind shear enhancing AC production is similar to 
that shown by Li and Schumacher [2011] over tropical 
Africa and South America. 

[ 35 ] RH profiles (Figure 8c) do not show much difference 
in the upper troposphere, but instead larger ACs are associ- 
ated with slightly lower RH in the midtroposphere between 
300 and 600 mb. One possible explanation is that for DCSs 
with very little AC, the NARR profiles in the AC region are 
more representative of the CC and SR region, and, therefore, 
they yield higher RH than systems with larger AC, which 
would be further away from active convection. Another 
possible reason is that errors and/or biases in upper tropo- 
spheric DCS moisture profiles produced by NARR prevent 
any significant signals to be obtained. 

[36] Vertical pressure velocity (Figure 8d) shows rising 
motion throughout the middle to upper troposphere for all 
systems investigated. Systems with smaller ACs (R A c <80 
km) tend to have stronger vertical motion, possibly for the 
same reason suspected for the greater RH, that the data are 
sampled closer to the convection. For systems with larger 
ACs (R A c >80 km) the rising motion between 1 50 to 600 mb 
increases with AC radius. We note that results from the 
vertical motion analysis should be treated with caution 
because comparisons between NARR, sounding and other 
reanalysis data sets show significant differences in vertical 
motion during the warm season at the Atmospheric Radia- 
tion Measurement (ARM) southern great plains site in 
Oklahoma [see Kennedy et al., 2011, Figure 4]. However, 
zonal and meridional wind fields from NARR compare 
favorably with ARM soundings, and RH also agrees well 
until above 300 mb where NARR has a 10-15% bias 
[Kennedy et al., 2011]. 

[ 37 ] Since the eventual size of the DCS is undeniably 
influenced by the lifetime of the system (Figure 4), the 
maximum size of AC attained is also related to its lifetime. 
Machado and Laurent [2004] found that the cloud area 
expansion rate at the initial stage of convection is a good 
indicator for system lifetime. The area expansion rate is 
defined as: 


where A is the convective system area defined by satellite 
T ir <235 K, V • V is the upper level wind divergence, and Q 
is the water content of the convective system [Machado and 
Laurent, 2004, equation (8)]. They suggested that a rapidly 
growing DCS is a result of strong liquid water condensation 
and vertical mass flux, along with increasing upper air wind 
divergence. 

[38] To test this hypothesis, area expansion rate and two 
radar-derived convective intensity indices of the initial hour 
are plotted as a function of total system lifetime in Figure 9. 
It is apparent that the area expansion rate exponentially 
increases with system lifetime until ~6 h, but the relation- 
ship no longer holds for lifetimes longer than 8 h. Rain rate 
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Figure 8. Composite NARR atmospheric state profiles as a function of the maximum AC radius attained 
during developing and mature stages of convection, (a) Zonal wind, (b) meridional wind, (c) relative 
humidity, and (d) vertical pressure velocity (negative is upward). The profiles are obtained during devel- 
oping and mature stages of each system and averaged within all anvil cloud regions; the numbers in the 
parentheses indicate the number of systems for that profile. 


and 30 dBZ echo top height show similar trends, except that 
the average 30 dBZ echo top height continues increasing 
with lifetime up to 10 h while the relationship with con- 
vective rain rate becomes ambiguous for lifetime longer than 
6 h. Nevertheless, these relationships support the findings of 
Machado and Laurent [2004] that rapid growth of convec- 
tive systems in the initial hour is a result of strong water 
condensation associated with vigorous convective forcing 
(e.g., strong low-level convergence, upper level divergence 
and vertical updraft). The results suggest that the lifetime of 
short to medium (less than 6 h) DCSs is strongly affected by 
the intensity of convection at the initial stage. However, for 
DCSs with longer lifetimes the vigor of initial convection no 
longer explains the system lifetime. Other factors such as a 
convectively favorable environment and the organization of 
convection must play a role to sustain DCSs lasting longer 
than 6 h. 

[ 39 ] To further investigate the impact of environmental 
factors on system lifetime, NARR reanalysis fields were 
composited based on the lifetime of systems as illustrated in 
Figure 10. NARR profiles within the CC and SR regions 
were averaged to obtain the composite. Compositing NARR 
within other regions of DCSs does not show significant 
differences, suggesting the signals are robust among differ- 
ent system lifetimes. Similar to Figure 8, NARR profiles 


during the developing and mature stages of the tracked 
DCSs are used for the composite. To maintain a robust 
classification of DCS life stage, only systems with lifetime 
longer than 3 h are considered. 



Figure 9. Average convective intensity at the initial hour 
as a function of the convective system lifetime. Area expan- 
sion rate ( A e ) is calculated using cloud area with T IR <235 K, 
and convective rain rate and 30 dBZ echo top height are 
derived from radar data. The shaded area indicates the num- 
ber of systems less than 1% of the total sample number. 
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Figure 10. Same as Figure 8, except composited as a function of system lifetimes. The profiles are 
obtained during developing and mature stages of each system and averaged within the convective and 
stratiform rain region. 


[ 40 ] While zonal wind profiles do not appear much dif- 
ferent between various system lifetimes, meridional wind 
shows stronger deep layer shear associated with longer- 
lasting systems. The bulk wind shear (between 150 and 
600 mb) values are similar between 4 and 6 h lifetime sys- 
tems (^20 m s _1 ), but monotonically increase to 27 m s -1 
for 9 h lifetime systems (~40% increase). RH differences 
between various system lifetimes are largest in the mid- 
troposphere. Similar to wind shear, RH show little differ- 
ences between 4 and 6 h lifetimes, but the mean RH 
value between 300 and 700 mb increases from 39% for 
6 h systems to 58% for 9 h systems (^50% increase). 
The differences between these profiles suggest that longer- 
lasting systems (more than 6 h) tend to occur in a convectively 
favorable environment (i.e., stronger wind shear and higher 
humidity), where SR region can fully develop several hours 
after the initial vigorous convection (Figure 6). While it is 
possible that high RH is partially caused by detrainment or 
evaporation/sublimation in the CC/SR region (i.e., NARR is 
contaminated by the convective signal), we argue that this is 
not the main reason for the difference because averaging 
the NARR profiles in the thin AC region (i.e., furthest away 
from the influence of convection) lead to similar conclusion. 

[ 41 ] Stronger deep layer wind shear and moisture associ- 
ated with longer-lived convective systems are consistent 


with previous studies of MCS environments in the United 
States [Cohen et al., 2007; Coniglio et al., 2007, 2010]. 
Vertical pressure velocity shows that short-lived systems are 
mostly associated with neutral to slightly downward motion 
throughout the troposphere, while longer-lived systems are 
related to upward motion in the middle to upper troposphere, 
consistent with mesoscale circulations in the SR region. The 
vertical velocity profile for 9 h lifetime systems shows 
transitioning from rising to sinking motion at around 500 mb. 
This is possibly due to cooling in the extensive SR region 
associated with long-lasting systems (Figure 5e) creating 
downward motion that eventually dominates the large-scale 
rising motion in the reanalysis. 

[ 42 ] Combining the evidence from Figures 6, 9, and 10, 
the following interpretation of DCS lifetime is provided. For 
short to medium systems (lifetimes less than 6 h), the life- 
time is mainly attributed to the intensity of the initial con- 
vection. More rapid expansion at the initiation stage signifies 
stronger initial convective forcing, resulting in longer system 
lifetime. For systems that last longer than 6 h without 
merging with other systems, vigorous initial convection 
alone is not sufficient to explain their lifetime. A moist 
midtroposphere with stronger deep layer wind shear is nec- 
essary to allow convection to fully deepen, providing an 
environment favorable for the growth of detrained water 
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condensates by deposition, releasing latent heat in the pro- 
cess and supporting further development of SR and AC 
region. Full development of the SR and AC regions takes a 
few hours longer after the initial intense convection, thus 
prolonging the lifetime of convective systems beyond 6 h. 
By incorporating radar and satellite data in a Lagrangian 
approach, and building on the findings of Machado and 
Laurent [2004] and Del Genio et al., [2012], this study 
provides interpretation of convective organization and the 
relationship between DCS lifetime, cloud properties, and 
their dynamic and thermodynamic environments. 

5. Summary and Conclusion 

[ 43 ] In this study, an automated satellite tracking method 
was used in conjunction with a multisensor hybrid classifi- 
cation to analyze the life cycle of precipitation and cloud 
structures associated with midlatitude deep convective sys- 
tems (DCSs) in a Lagrangian framework during two warm 
seasons (May- August 2010-2011) over the central United 
States. 

[ 44 ] On average, the frequency of occurrence for all 
DCSs tracked by the automated satellite tracking method 
is -'-43 day' 1 during the study period. The average DCS 
cloud fraction accounts for 25% of total clouds as observed 
by GOES satellite, with 8% of total clouds from precipitating 
deep clouds. The majority of the tracked systems have life- 
times shorter than 14 h, and the analysis is confined within 
the southern and central United States where these sys- 
tems are well covered by both NEXRAD radars and 
GOES satellite. The number of systems decreases expo- 
nentially with increasing lifetime, and ~51% of the sys- 
tems have a lifetime of 3 h or longer. Anvil clouds (AC) 
dominate the DCS area, covering an area about three 
times that of stratiform rain (SR) area and an order of 
magnitude that of the convective core (CC). SR accounts 
for 72% of the precipitation area, similar to the tropics- 
wide average SR area fraction reported in other studies. 
Maximum system size attained during their life cycles 
correlates with system lifetime. Merged and nonmerged 
systems mainly differ in lifetime and maximum size, with 
majority of the nonmerged systems lasting shorter than 
10 h. 

[ 45 ] Statistical analysis of the DCS evolution shows that 
longer-lived DCSs are consistently larger with more exten- 
sive stratiform rain regions and anvil cloud shields. Com- 
posites of the maximum SR and AC area lag peak 
convective intensity (Cl), where the lag increases linearly 
with system lifetime, and varies from ~1 h for short-lived 
systems to more than 3 h for long-lived ones. For short to 
medium systems (lifetimes less than 6 h), the lifetime is 
mainly attributed to the intensity of the initial convection. 
For systems that last longer than 6 h without merging with 
other systems, vigorous initial convection alone is not suf- 
ficient to explain their lifetime. These longer-lasting systems 
are associated with up to 50% higher midtropospheric rela- 
tive humidity and up to 40% stronger middle to upper tro- 
pospheric wind shear. Such environments allow continuous 
growth of detrained hydrometeors by deposition after the 
initial intense convection, releasing latent heat in the process 
and supporting further development of SR and AC region. 
The increased staggered timing with lifetime between 


maximum SR/AC area and peak convective intensity indi- 
cates that changes in the overall diabatic heating structure 
associated with the transition from CC to SR and AC could 
help maintain the systems through generating mesoscale 
ascent, prolonging the system lifetime beyond 6 h by sus- 
taining SR and AC development [Del Genio et al., 2012], 

[46] Factors that are most relevant to the convective AC 
production were investigated and results show that the areal 
coverage of thick AC is strongly correlated with the size of 
CC, updraft strength, and the SR area. After the initial 
intense convection in the developing stage, AC from 
detrainment of CC or remnants of SR precipitation continues 
spreading out with upper level wind. Stronger upper tropo- 
spheric wind speed and middle to upper tropospheric wind 
shear contribute directly to producing large areas of AC until 
systems mature. For example, the average wind speed (150 
~ 200 mb) and wind shear (150 to 600 mb) for systems with 
AC radius >120 km are 24% and 20% higher, respectively, 
than those with AC radius = 20 km. When convection ceases 
and moisture supply to the AC region is cut off at the 
decaying stage, AC dissipates within a few hours. The 
influence of moisture on AC growth appears not directly 
through higher RH in the upper troposphere that favors 
growth of detrained ice, but rather indirectly through higher 
midtropospheric moisture, along with stronger deep layer 
shear that favor longer-lived and more intense DCSs, to 
produce larger size of convective AC. 

[ 47 ] Analyses performed in this study provide a frame- 
work for comparison with long-term large-scale cloud 
resolving model (CRM) simulations to systematically eval- 
uate their performance in simulating convections at various 
scales. Classification and automated cloud tracking can be 
applied on proxy observation data set from model output 
fields (e.g., simulated radar reflectivity and cloud top 
brightness temperature), allowing comparison of precipita- 
tion and cloud structure statistics directly with observations 
in the Lagrangian framework. The evolution of the CC/SR/ 
AC, their lifetime and associated environments, and factors 
that impacts AC production found in this study should also 
be compared with more CRM simulated cases to diagnose 
the mechanisms responsible for transition from CC to SR 
and AC, and provide guidance for parameterizing these 
processes in Global Climate Models. A recent field cam- 
paign conducted in Oklahoma (Midlatitude Continental 
Convective Clouds Experiment, or MC 3 E, http://campaign. 
ann.gov/mc3e/ and http://mc3e.nsstc.nasa.gov/) incorpo- 
rated high-resolution multiwavelength radar systems, sur- 
face and in situ aircraft measurements, and intensive 
atmospheric sounding observations. Data sets generated 
during MC 3 E have great potential for complementing the 
large-scale DCS morphology study performed in this paper 
through comprehensive comparison of the microphysical 
processes with CRM simulations. 
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